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Abstract 

New CP 1 -soliton behaviour on a flat torus is reported. Defined by 
the Weierstrass elliptic function and numerically-evolved from rest, 
each soliton splits up in two lumps which eventually reunite, divide 
and get back together again, etc.. This result opens up the question 
of fractional topological charge. 

1 Introduction 

The CP 1 model in (2+1) dimensions appears as a low dimensional 
analogue of non-abelian gauge field theories in four dimensional space- 
time. This analogy relies on common properties like conformal invari- 
ance, existence of topological solitons, hidden symmetry and asymp- 
totic freedom. Amongst various applications, CP 1 models have been 
used in the study of the quantum Hall effect and high-T c superconduc- 
tivity. In differential geometry, the soliton-solutions of CP 1 models 
are known as harmonic maps, a rich industry of research on its own. 

The classical (2+0)-dimensional CP 1 or non-linear 0(3) model on 
the extended plane 3ft 2 U {oo} « 52, where the soliton solutions are 
harmonic maps £2 |— > S%, has been amply discussed in the literature 
0) @|- I n (2+1) dimensions the model is not integrable, and the study 
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of its dynamics is done with the aid of numerical simulations. Due to 
the conformal invariance of the theory on the plane, the 0(3) solitons 
are unstable in the sense that they change their size under any small 
perturbation, either explicit or introduced by the discretisation pro- 
cedure. It can make the solitons shrink indefinitely and, when their 
width is comparable to the lattice spacing, the numerical code breaks 
down ||. However, such instability can be cured by the addition of 
two extra terms to the lagragian [|j . The first one resembles the term 
introduced by Skyrme in his nuclear model in four dimensional space- 
time ||, and the second one is a potential term. The fields of the 
planar Skyrme model (skyrmions) produce stable lumps which repel 
each other when started off from rest j|, ||]. 

In a recent paper we considered both the pure and modified 
CP 1 schemes imposing periodic boundary conditions, which amounts 
to defining the system on a torus T^. The corresponding soliton con- 
figurations are harmonic maps T2 1— > S2. In ]7|] we found (using the 
Weierstrass' cr(z) function to define the solitons) that in contradis- 
tinction with the familiar theory on S2, the toroidal model: • has no 
analytical single-soliton solution [this is because elliptic functions, in 
terms of which the toroidal solitons must be expressed, are at least of 
the second order; or, in the language of differential geometry, because 
genus(torus)=l]; • needs only a Skyrme term to stabilise the solitons 
(thus the lagrangian retains its 0(3) invariance: on S2, the latter is 
broken by the additional potential term); • does not require a damp- 
ing set-up for the numerical simulation (a radiation-absorbing device 
is implemented for the model evolved on the compactified plane in 
order to prevent the reflection of kinetic waves from the boundaries); 
• has perfectly static skyrmions when their initial speed vo is zero 
(as already pointed out, on £2 they move away from each other for 
Vq = 0); • possesses no critical velocity below which the skyrmions 
scatter back-to-back in head-on collisions. They always scatter at 
right angles provided vq 7^ 0. Also, on T2 the skyrmions scatter any 
number of times (multi-scattering), as they keep encountering each 
other in the periodic grid. 

In the present work we continue the study of periodic CP 1 con- 
figurations, limiting ourselves to those with vq=0 in the topological 
charge- two sector. Defining the solitons through the elliptic func- 
tion p(z) of Weierstrass we will see new soliton behaviour, where each 
skyrmion-lump divides itself in two smaller components that glue back 



2 



together, split and reunite again and so on. 



2 Periodic skyrmion model 

Our model is given by the lagrangian density 

\dtWl-mwl \d z W\ 2 a 2 

£ - (TTWW 2 jr+ww il 1 " 11 } ' ( } 

z = x + iy e T 2 , which is the pure CP 1 model plus an additional 
Skyrme, #i-term {0\ e The complex field W obeys the periodic 

boundary condition 

W[z + (m + in)L] = W(z), Vi 5 (2) 

where m, n = 0, 1, 2, ... and L is the size of a square torus. The static 
solitons (skyrmions) are elliptic functions which may be written as 

W = A p(z — a) + b, \,a,beZ, (3) 

p{z) being the elliptic function of Weierstrass. Within a fundamental 
cell of length L, p possesses the expansion || 

p(z) = z- 2 + i 2 z 2 + £ 3 z 4 + + i 3 z 2 i~ 2 + £j e ». (4) 

This function is of the second order, hence (|3|) represents solitons 
of topological index 2. Note that (|3|) is an approximate solution of 
the model ([!]), except in the pure CP 1 limit (#i=0) where it exactly 
solves the corresponding static field equation. Therefore, we expect 
our solitons to evolve only for a non-zero Skyrme parameter. 
In reference we computed the periodic solitons through 

employing a subroutine that numerically calculates cr(z). Via the for- 
mula below ||, in this paper we use the same subroutine to compute 
P(z) ■ 

d 2 

p{z) = -^2 M a 0)]> ( 6 ) 

where the Laurent expansion for a reads 

00 

a(z) = Y / c j z Aj+ \ Cj e5R. (7) 

j=0 
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3 Basic numerical procedure 



We treat configurations of the form (||) as the initial conditions for 
our time evolution, studied numerically. Our simulations run in the 
^-formulation of the model, whose field equation follows from the la- 
grangian density (|l]) with the help of the stereographic projection 

W = ^±^, (8) 

1-03 

where the real scalar field <j> = (<f>i, fo, ^3) satisfies 4>.(f) = 1. 

We compute the series (0) up to the fifth term, the coefficients cj 
being in our case negligibly small for j > 6. We employ the fourth- 
order Runge-Kutta method and approximate the spatial derivatives by 
finite differences. The laplacian is evaluated using the standard nine- 
point formula and, to further check our results, a 13-point recipe is also 
utilised. Our results showed unsensitiveness to either method. The 
discrete model evolves on a 200 x 200 periodic lattice (n x = n y = 200) 
with spatial and time steps 5x=5y=0.02 and <5i=0.005, respectively. 
The size of our fundamental, toroidal network is L = n x x Sx = 4. 

Unavoidable round-off errors gradually shift the fields away from 

the constraint (f>.<j) = 1. So we rescale <j> — > <j)/ J cf.cf every few itera- 
tions. Each time, just before the rescaling operation, we evaluate the 
quantity [i = — 1 at each lattice point. Treating the maximum 
of the absolute value of /i as a measure of the numerical errors, we 
find that max|/x I w 10" 8 . This magnitude is useful as a guide to de- 
termine how reliable a given numerical result is. Usage of an unsound 
numerical procedure in the Runge-Kutta evolution shows itself as a 
rapid growth of max|/i|; this also occurs, for instance, in the 0(3) limit 
(#1 = 0) when the unstable lumps of energy become infinitely spiky. 



4 Splitting lumps 

The energy density associated with our degree-2 solitons 
W = Xp{z — a) + b can be read-off from the lagrangian (|l|). The 
pole a determines the position of the energy humps on the basic cell; 
we shall take a=(2.025, 2.05). The quantity b defines the distance be- 
tween the solitons; for 6=0 we have W = p{z—a) (use A=l throughout, 
for simplicity), whose energy density gives indistinguishable lumps on 



4 



top of each other, as depicted in figure [j] (top-left). The value b=l 
gives two lumps separated along the ordinates [figure [j] (top-right)], 
whereas b=-l positions them along the abscissas. A pure imaginary 
b places our extended structures on a diagonal bisecting the toroidal 
grid [for a=(2,2) the bisection is exact, but this is not a good numer- 
ical value because z, with spatial steps of 0.02, blows W at z = a], 
and b with non-zero real and imaginary parts situates the lumps in an 
arbitrary diagonal of the cell. These set-ups are true regardless of 9i, 
which we have put equal to 0.001. 

Our numerical simulations show that the skyrmions are stable. 
Their stability is reflected on the left-hand side of the nether half of 
figure |l[ which shows the evolution of the the maximum value of the 
system's total energy density (E max ) for 6=0,1. In the 0(3) limit the 
lumps are no longer stable, as can be appreciated from the bottom- 
right graph of figure [l|. In this expected, the solitons remain 
static with the passing of time. 

But in the stable, Skyrme situation, the lumps evolve in novel fash- 
ion. Let us first consider the configuration when the extended entities 
are on top of each other at t = 0. As time elapses, the system splits 
in four equal lumps, each progressing towards its nearest lattice cor- 
ner. There they meet and coalesce, for all corners are nothing but the 
same point. Then the system splits up once more and the 'fractional- 
skyrmions' make their way back to the centre of the nett, in a cycle 
that repeats itself indefinitely. The foregoing event is illustrated in the 
superior half of figure [2|, with the trajectory of the four energy peaks in 
the x — y plane. The accompanying 3-D picture captures the moment 
when the skyrmion quartet, having concurred at the corners and split 
afresh, begin to motion towards the centre of the network. Worthy of 
remark is that the trajectory in question resembles the usual head-on 
collision course and subsequent 90° scattering of two solitons, in spite 
of ours being energy chunks with no initial velocity. Also, we have 
thoroughly verified that the topological number is 2 all along our nu- 
merical evolution, which suggests that each 'fractional soliton' carries 
a topological charge of 1/2. 

Two CP^solitons on top of each other on S2 have been numeri- 
cally studied in Q. Such initially coalescing objects (in the Skyrme 
variation of the model) were found to move away from one another, 
two evolving lumps in mutual repulsion. On T2, however, we have 
observed a qualitatively different phenomenon. 
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A more involved trajectory occurs for two initially well-separated 
skyrmions. In the bottom-left plot of figure the labels a — g indicate 
the itinerary of one of the entities (call the corresponding symmetrical 
points a! — g'). At t = a full lump is at point a but it soon halves 
under the numerical simulation. One of its fractional offspring moves 
following curve b, whereas its counterpart proceeds in the opposite 
sense. At x = = 4 they get back together into one full structure, 
which runs vertically up before separating anew. One of these com- 
ponents cruises along c and, at site d, reunites back with its peer 
travelling from the left. Before dividing itself according to curve e, 
the skyrmion is seen to shift towards the centre, as one can tell from 
the small leg connecting curves d-e (the full lump started at a' under- 
goes a similar process). The bottom-right diagram of figure || exhibits 
the skyrmions heading centrewards from d and d'. Thence our system 
continues through f — g — h and returns to its t = coordinates. Ob- 
serve that at g a half-soliton from lump a undistinguishably coalesces 
with a half-soliton from lump a' , so actually we do not know which bit 
is ascending (descending) along h (h'). We terminated our simulations 
when a like cycle was about to commence, as evidenced by the small 
vertical lines emerging from a, a'. 

We underline that our research has been constraint to systems 
with zero inital speed. Important mathematical aspects of CP 1 soli- 
tons given by equation (^) [the 0(3) case only] have been recently 
analysed in ||], where the dynamics of the solitons was carried out 
using the geodesic approximation. The presence of four energy peaks 
rather than two is therein discussed as well. Through our numerical 
approach, we are planning to study ourselves the dynamics of soliton 
configurations defined by (g). 

Finally, note that the study carried out in this article shows that 
both repulsive and attractive interactions are operating between the 
skyrmions. 
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5 Concluding remarks 



The CP 1 model in (2+1) dimensions is variegated. More so its stable, 
skyrmionic version on T2, which in this work has been shown to possess 
qualitatively different features as compared to the familiar model on 
the compactified plane. 

A particularity of our periodic solitons is that they have no ana- 
lytic representative of degree one, limitation dictated by their elliptic 
nature. In (2+0) dimensions, the CP 1 model on S2 is known to have 
soliton solutions in all topological classes. 

Another peculiarity of the toroidal model, one herein discovered, 
is that the properties of the skyrmions depend on the elliptic function 
used to define them. Thus, skyrmion fields expressed in terms of 
c(z), equation (||), evolve differently than those expressed through 
p(z), equation (||). In the former case, for vq = 0, the associated 
chunks of energy stay still in their initial positions as time goes by. 
In the latter case, the system splits up in four lumps that stroll the 
network, comeback together, etc., acted upon by repulsive-attrative 
forces. Consequently, it would be interesting to analyse the solitons 
on T2 by employing alternative elliptic functions. Although in another 
context, a step in this direction has already been taken in |K|, where 
instantons on a torus are displayed in terms of a Jacobian elliptic 
function. Let us remind that in S% no new traits arise from casting two- 
soliton configurations in different ways, e.g., W = z 2 , z~ 2 , ^Z.^z-d) . 

That the topological charge of our splitting system is very well con- 
served (=2) as time progresses, invites speculation on whether each 
'fractional-soliton' carries a non-integer degree. Clearly, further re- 
search on this matter is required. 

Investigation on the appealing question of collisions between the 
solitons (H) is currently under way. Note that already on T2, the 
fields (|5|) have yielded -in the Skyrme case- the outcome of always 
scattering off at 90° when impinged with a non-zero initial speed Q. 
On S2, unsimilarly, the existence of a critical speed, below which the 
skyrmions scatter at 180° to the initial direction of motion, has long 
been a landmark of the planar model . 
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Initial configuration for b=0 



Initial configuration for b=1 




20 40 60 5 10 15 

t t 

Figure 1: Energy density configurations at t — and the evolution of their 
peaks. The bottom-right graph corresponds to #i=0, when the lumps are 
unstable and shrink non-stoppingly. 
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Figure 2: Above: Trajectory of the skyrmions initially on top of each other. 
They split up in four lumps heading to the corners where they coalesce and 
break-off again, moving back to the centre of the lattice, as in the illustration 
for £=51. Below: The initially separated skyrmions also divide each in two, 
but transit more complicated paths; the signs a — h refer to one of the 'half- 
lumps'. The t=30 picture is shortly after the fractional progeny have reunited 
at d (and at its symmetrical point) and begun to travel centrewards. 



9 



Acknowledgements 



I thank J. M. Speight for drawing my attention to p-defined 
solitons, and W. J. Zakrzewski for enlightening conversations. 
The financial support of Universidad del Zulia is acknowledged. 

References 

[1] Eichenherr H. (1976) Nucl. Phys. B146 215 

[2] Perelomov A. M. (1981) Physica 4D 1 

[3] Leese R. A., Peyrard M. and Zakrzewski W. J. (1990) Non- 
linearity 3 387 

[4] Leese R. A., Peyrard M. and Zakrzewski W. J. (1990) Non- 
linearity 3 773 

[5] Skyrme T. H. R. (1962) Nucl. Phys. 31 556 

[6] Sutcliffe P. M. (1991) Nonlinearity 4 1109 

[7] Cova R. J and Zakrzewski W. J. (1997) Nonlinearity 10 1305 

[8] Goursat E. (1916) Functions of a complex variable Dover 
Publications. Alternatively, see: Erdelyi A. et. al. (1953) 
Higher transcendental functions vol II Mc Graw Hill 

[9] Speight J. M. Lump dynamics in the CP 1 model on the 
torus, University of Texas at Austin preprint, Dept. of Math, 
Austin, Texas, 78712, U.S.A 

[10] Sutcliffe P. M. (1995) Nonlinearity 8, 411 



10 



